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ABSTRACT 

We present deep, high-resolution imaging of the nearby Fanaroff- Riley Class I (FRI) 
radio galaxies NGC 193, B2 0206+35, B2 0755+37 and M84 at frequencies of 4.9 and 
1.4 GHz using new and archival multi-configuration observations from the Very Large 
Array. In addition, we describe lower-resolution observations of B2 0326+39 and a re- 
analysis of our published images of 3C 296. All of these radio galaxies show twin jets 
and well-defined lobes or bridges of emission, and we examine the common properties 
of this class of source. We show detailed images of total intensity, brightness gradient, 
spectral index, degree of polarization and projected magnetic-field direction. The jet 
bases are very similar to those in tailed twin-jet sources and show the characteristics 
of decelerating, relativistic flows. Except on one side of M 84, we find that the jets 
can be traced at least as far as the ends of the lobes, where they often form struc- 
tures which we call "caps" with sharp outer brightness gradients. Continuing, but less 
well collimated flows back into the lobes from the caps can often be identified by 
their relatively flat spectral indices. The lobes in these radio galaxies are similar in 
morphology, spectral-index distribution and magnetic-field structure to those in more 
powerful (FRII) sources, but lack hot-spots or other evidence for strong shocks at the 
ends of the jets. M 84 may be an intermediate case between lobed and tailed sources, 
in which one jet does not reach the end of its lobe, but disrupts to form a "bubble". 
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1 INTRODUCTION 

Relativistic jets are the primary channel of energy loss from 
accreting supermassive black holes in many radio galaxies. 
They also have a major impact on their surroundings and 
act as accelerators of the most energetic photons (and per- 
haps hadrons) we observe. The present paper forms part of a 
study of jet physics in nearby, low- luminosity radio galaxies, 
spe cifically those with FRI morphology (Fanar off &: Rilevl 
1 1974 ). We have developed a sophisticated model of FRI jets 
as relativistic, symmetrical, axisymmetric flows. By fitting to 
deep, high-resolution radio images in total intensity and lin- 
ear polarization, we have determined the three-dimensional 
variations of v elocity, emissivity and magnetic- field ordering 
in five sources (|Laing fc Bridle. 2002a: Canvin fc Laing:.2004 : 



ICanvin et al.ll2005l : iLaing et al.ll2006bh . We have shown that 
FRI jets decelerate from relativistic (/3 — v/c ~ 0.8) to sub- 
relativistic speeds on scales of a few kpc and that they are 
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faster on-axis than at their edges, as expected if they entrain 
external material. 



The physics of boundary-layer entrainment must de- 
pend on the composition and density of the surround- 
ing medium, and in particular on whether the jets propa- 
gate in direct contact with the intergalactic medium (IGM) 
or are surrounded by lobes consisting primarily of tenu- 
ous and at least partially relativistic plasma. Of the five 
sources we have modelled, three have plumed or tailed 
outer structures wherein most of the extended emission ap- 
pears to lie further from the active nucleus than t he nar - 
rower jets: 3 C31 (iLaing k B ridle 200 2al: iLaing et al.ll20Q8h . 
B2 1553+2 4 (ICanvin fc Laing. .2004: .Young et al.ll2005l ) and 
NGC 315 (ICanvin et al.ll2005l : iLaing et al.ll2006'al ). We pre- 
sume that their jets are i n direct contact wit h the IGM. On 
the other hand, 3C 296 (iLaing et all l2006bl ) has two lobes 
with well-defined outer boundaries and a diffuse bridge of 
emission around the jets (at least in projection ). A fur- 
ther source, B2 0326+39 (|Canvin fc Laingj l2004l ). clearly 
has lobes, but it was unclear from published observations 
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(|Bridle et al.lll99ll ) whether these lobes surround the inner 
jets. 

Our be st-fit mod el for the jets in the lobed FRI source 
3C 296 (Lai ng et alJ [2006b) is unusual in that it shows a 
very large transverse velocity gradient across the jet except 
very close to the nucleus, the ratio of edge to central ve- 
locity in this jet being ^0.1, compared with values ^0.4 
for the jets in three tailed sources (the value for 0326+39 
is poorly determined) . It is therefore of interest to examine 
whether jets in other FRI sources whose lobes entirely sur- 
round them appear to decelerate differently with distance 
from the nucleus than those in tailed FRI sources, or those 
in sources whose lobes may not extend all the way back 
to the nucleus, leaving the inner jets unshielded. Projection 
complicates interpretation of individual sources (the lobes 
may appear superimposed on the jets even if they are not 
in physical contact) and it is not always straightforward to 
separate jet and lobe emission, but the differences between 
3C 296 and the rest are large. Modelling of the jets in a small 
number of lobed sources (which form the maj ority of com- 
plet e samples of low- luminosity radio galaxies; Parm a et al.l 
1 1996 ) should therefore be enough to decide whether there 
are systematic differences between the two classes. 

The primary aim of this paper is to present high-quality 
radio imaging of four lobed FRI sources whose jets are suit- 
able for modelling by our methods. The models themselves 
will be presented elsewhere. Our approach to jet modelling 
requires fitting to high-fidelity, deep images with linear res- 
olution <0.25kpc derived from mult i- configuration obser- 
vations at 4.9 GHz (C-band) or 8.5 GHz (X-band) with the 
Very Large Array (VLA) at the National Radio Astronomy 
ObservatorjQ. In order to correct the linear polarization for 
the effects of Faraday rotation (small at these frequencies), 
we also need to observe at several frequencies in the range 1.3 
- 1.7 GHz (L-band). High-quality images of depolarization, 
rotation measure and spectral index are useful byproducts. 
In this paper, we describe: 

(i) details of the observations, 

(ii) the source morphologies in total intensity at a range 
of resolutions, 

(iii) images of the spectral-index distributions and 

(iv) images of the degree of polarization and the apparent 
magnetic field direction, corrected for Faraday rotation. 

We also include high-resolution MERLIP^fl imaging for two 
of the sources. Fa raday rotation a nd de polarization are anal- 
ysed in detail by iGuidetti et al.l (|201lL and in preparation) . 



We also present an analysis o f the spectrum of 3 C 296 
which improves on that given by iLaing et al.l (|2006b|) and 
imaging of the large-scale structure of B2 0326+39 to trace 
its lobe emission closer to the nucleus than in earlier stud- 
ies. These data complete the documentation of the large- 
scale structures and spectral-index distributions for all of 
the lobed FRI sources whose jets we can currently model. 



Table 1. The sources. Col. 1: name (as used in this paper). 
Col. 2: alternative names. Col. 3: redshift. Col. 4: linear scale 
(kpcarcsec"-*^) for our adopted cosmology. Col. 5: reference for 
redshift. Col. 6: reference for earlier observations of large-scale ra- 
dio structure. 



Name 



Alternative 



kpc 



Ref 
O R 



NGC193 PKS 0036+03 0.0147 0.300 



0206+35 

0755+37 
M84 

3C296 
0326+39 



UGC408 
UGC1651 
4C 35.03 
NGC 2484 
3C 272.1 
NGC 4374 
NGC 5532 



0.0377 0.748 

0.0428 0.845 

0.0035 0.073 

0.0247 0.498 

0.0243 0.490 



9 

3 

9 
10 

7 
7 



Re ferences 
3 iFalco et al 



1 iBondi et al.l (l2000l). 2 iBridle et al.l (Il99lh . 
" I999l ). 4 iGiacintucci et al.l (l201lh . 5 



Bridk 



6 iLaing et al.l (l2006bl)~7 iMiller et al.l 

Morganti et al.l (Il987l ). 9 lOgando et al.l (|2008l ). 10 



[Laing 

120021 

iTrager et al.l (|2000l ) 



Section [2] describes the new observations and data re- 
duction. Section [3] presents our results for the sources indi- 
vidually, and Section |4] outlines the phenomenology of their 
jets and larger-scale emission as a prelude to modelling. Sec- 
tion [5] is a brief summary. 

We adopt a concordance cosmology with Hubble con- 
stant, Hq = 70kms~^ Mpc"\ Qa = 0.7 and Qm = 0.3. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The sources 

We selected four bright FRI sources for which available 
data suggested: (a) that full synthesis observations with the 
VLA would achieve signal-to-noise sufficient to image the 
linearly polarized emission from their counter-jets with sev- 
eral beamwidths resolution transverse to the radio features, 
as required by our modelling methods, and (b) that the jets 
have formed lobe-like structures rather than diffuse outer 
plumes. Three of the sources, NGC 193, B2 0206+35 and 
B2 0755+37 are analogous to 3C 296 in having well-defined 
outer boundaries. The fourth, M 84, is in some respects an 
intermediate case between the two classes, as we discuss be- 
low. We analysed a combination of new and archival datasets 
chosen to give good spatial-frequency coverage at two or 
three frequencies. 

In addition, we improv ed our low-resolution images of 
3C296 ([Laing et al.l l2006bl V Finally, we analysed shorter, 
low-resolution, archival VLA observations for B2 0326+39. 

Alternative names, redshifts, linear scales and references 
for all of the sources are given in Table [1] (we drop the B2 
from source names from now on). A journal of observations 
is given in Table [2n. 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

^ MERLIN is a UK National Facility operated by the University 
of Manchester at Jodrell Bank Observatory on behalf of STFC. 



^ iLaing et al.ll2006bl give a full description of the observations of 
3C 296; this is not repeated here. 
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Table 2. Journal of VLA observations. Col. 1: source name. Col. 2: VLA configuration. Col. 3: Date of observation. Col. 4: centre fre- 
quencies for the one or two channels observed (MHz). Col. 5 bandwidth (MHz). Col. 6: the on-source integration time scaled to an array 
with all 27 antennas operational. Col. 7: VLA proposal code. 



Source Config- 


Date 


u 




Aiy 


t 


Proposal 


uration 




[MHz] 




[MHz] 


[min] 


code 


NGC 193 A 


2007 Jun 02 


4885.1, 


4835.1 


50 


280 


AL693 


A 


2007 Jun 28 


4885.1, 


4835.1 


50 


90 


AL693 


A 


2007 Aug 22 


4885.1, 


4835.1 


50 


41 


AL693 


A 


2007 Aug 23 


4885.1, 


4835.1 


50 


88 


AL693 


B 


2007 Nov 05 


4885.1, 


4835.1 


50 


318 


AL693 


B 


2007 Nov 16 


4885.1, 


4835.1 


50 


121 


AL693 


C 


2008 May 24 


4885.1, 


4835.1 


50 


223 


AL693 


D 


2007 Mar 11 


4885.1, 


4835.1 


50 


53 


AL693 


A 


2007 Jun 28 


1365.0 




25 


83 


AL693 


A 


2007 Aug 22 


1365.0 




25 


39 


AL693 


A 


2007 Aug 23 


1365.0 




25 


97 


AL693 


B 


2007 Nov 16 


1365.0 




25 


148 


AL693 


C 


2008 May 24 


1365.0 




25 


61 


AL693 


0206+35 A 


2008 Oct 13 


4885.1, 


4835.1 


50 


486 


AL797 


A 


2008 Oct 18 


4885.1, 


4835.1 


50 


401 


AL797 


B 


2003 Nov 17 


4885.1, 


4835.1 


50 


254 


AL604 


C 


2004 Mar 20 


4885.1, 


4835.1 


50 


88 


AL604 


A 


2004 Oct 24 


1385.1, 


1464.9 


25 


189 


AL604 


B 


2003 Nov 17 


1385.1, 


1464.9 


25 


110 


AL604 


0755+37 A 


2008 Oct 05 


4885.1, 


4835.1 


50 


477 


AL797 


A 


2008 Oct 06 


4885.1, 


4835.1 


50 


383 


AL797 


B 


2003 Nov 15 


4885.1, 


4835.1 


50 


332 


AL604 


B 


2003 Nov 30 


4885.1, 


4835.1 


50 


169 


AL604 


C 


2004 Mar 20 


4885.1, 


4835.1 


50 


125 


AL604 


D 


1992 Aug 2 


4885.1, 


4835.1 


50 


55 


AM364 


A 


2004 Oct 25 


1385.1, 


1464.9 


12.5 


450 


AL604 


B 


2003 Nov 30 


1385.1, 


1464.9 


12.5 


160 


AL604 


C 


2004 Mar 20 


1385.1, 


1464.9 


12.5 


21 


AL604 


M84 A 


1980 Nov 09 


4885.1 




50 


223 


AL020 


A 


1988 Nov 23 


4885.1, 


4835.1 


50 


405 


AW228 


A 


2000 Nov 18 


4885.1, 


4835.1 


50 


565 


AW530 


B 


1981 Jun 25 


4885.1 




50 


156 


AL020 


C 


1981 Nov 17 


4885.1 




50 


286 


AL020 


C 


2000 Jun 04 


4885.1, 


4835.1 


50 


138 


AW530 


A 


1980 Nov 09 


1413.0 




25 


86 


AL020 


B 


1981 Jun 25 


1413.0 




25 


29 


AL020 


B 


2000 Feb 09 


1385.1, 


1464.9 


50 


30 


AR402 


0326+39 D 


1997 Dec 13 


4885.1, 


4835.1 


50 


11 


AR386 


D 


1997 Dec 16 


4885.1, 


4835.1 


50 


32 


AR386 


C 


1998 Dec 4 


1464.9, 


1414.9 


50 


11 


AR386 


D 


1997 Dec 13 


1464.9, 


1385.1 


50 


6 


AR386 


D 


1997 Dec 16 


1464.9, 


1385.1 


50 


9 


AR386 



2.2 VLA data reduction 

The VLA data listed in Table [2] were calibrated and imaged 
using the AiPS software package, following standard proce- 
dures with a few additions. The flux-density scale was set us- 
ing observations of 30 286 or 30 48 and (except for 0326+39) 
the zero-point of E- vector position angle was determined us- 
ing 30 286 or 30 138, after calibration of the instrumental 
leakage terms. The main deviations from standard methods 
were as follows. 

Firstly, we used the routine blcal to compute clo- 
sure corrections for the 4. 9- GHz observations. This was re- 
quired to correct for large closure errors on baselines be- 
tween E VLA and VLA an tennas in observations from 2007 
onwards (|VLA StafJl201Ql ). but also improved a number of 



the earlier datasets. Whenever possible, we included obser- 
vations of the bright, unresolved calibrator 30 84 for this 
purpose; if it was not accessible during a particular observ- 
ing run, we used 30 286. We found that it was not adequate 
to use the standard calibration (which averages over scans) 
to compute the baseline corrections, as phase jumps during 
a calibrator scan caused serious errors in the derived correc- 
tions. We therefore self- calibrated the observations in ampli- 
tude and phase with a solution interval of 10 s before running 
BLCAL. We assumed a point-source model for 30 84 and the 
well-determined clean model supplied with the AiPS distri- 
bution for 30 286. 

Secondly, we imaged in multiple facets to cover the inner 
part of the primary beam at L-band and to image confusing 
sources at large distances from the phase centre in all bands. 
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Table 3. Resolutions and noise levels for the images used in this paper. Col. 1: source name. Col. 2: resolution (FWHM, in arcsec). 
Col. 3: frequency, in MHz. Col. 4: VLA configurations used in the image (not relevant for the two MERLIN observations). Cols 5 and 
6: deconvolution method for / and Q/U images, respectively (MR: multi-resolution CLEAN; CL: single-resolution CLEAN; ME: maximum 
entropy). Cols 7 and 8: off-source noise levels, crj is the noise level on the / image; ap the average of the noise levels for Q and U . All 
noise levels were determined before correction for the primary beam response and larger v alues are appropriate at significant distances 
from the field centre. Col. 9: the approximate maximum scale of structure imaged reliably (|Ulvestad. Perlev fc Chandlej 120091 ). Cols 10 
and 11: fiux densities measured from the image and derived from single-dish measurements, respectively. The singl e-dish fiux densities 
were taken from the references in Col. 12. They have been corrected as necessary to put them on the standard scale of lBaars et al.l (|l977l ') 
and interpolated to our observing frequencies. 



Source 


FWHM 


Freq 


Config- 


Method 


rms noise level 


Max 




/int/Jy 


Reference 




[arcsec] 


[MHz] 


urations 


/ QU 


[/iJybeam"-*^] 
aj ap 


scale 

[arcsec] 


Image 


SD 





NGC 193 



0206+35 



0755+37 



M84 



3C296 



0326+39 



0.45 

1.35 

1.60 

1.60 

4.05 

4.05 

0.16 

0.35 

1.20 

1.20 

1.20 

1.20 

4.50 

4.50 

0.14 

0.40 

1.30 

1.30 

1.30 

1.30 

4.00 

4.00 

4.00 

4.00 

0.40 

1.65 

1.65 

4.5 

4.5 

4.5 

4.5 

5.5 

5.5 

18.0 

18.0 



4860.1 

4860.1 

4860.1 

1365.0 

4860.1 

1365.0 

1658.0 

4860.1 

4860.1 

1464.9 

1385.1 

1425.0 

4860.1 

1425.0 

1658.0 

4860.1 

4860.1 

1464.9 

1385.1 

1425.0 

4860.1 

1464.9 

1385.1 

1425.0 

4860.1^ 

4860." 

1413.0 

4860.1^ 

1464.9 

1413.0 

1385.1 

8460.1 

1479.0 

4860.1 

1425.0 



ABCD 

ABCD 

ABCD 

ABC 

ABCD 

ABC 

MERLIN 

ABC 

BC 

AB 

AB 

AB 

BC 

AB 

MERLIN 

ABCD 

BCD 

ABC 

ABC 

ABC 

BCD 

ABC 

ABC 

ABC 

ABC 

ABC 

AB 

C 

B 

AB 

B 

ABCD 

BCD 

D 

CD 



ME 

MR 

MR 

MR 

CL 

MR 

CL 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

CL 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

MR 

CL 

CL 



CL 

CL 

CL 

MR 

CL 

MR 

MR 
MR 
MR 
MR 



MR 
MR 
MR 
MR 

MR 
MR 
MR 
MR 
MR 
MR 
MR 
MR 
MR 
MR 
MR 



8.6 
7.1 
7.5 
36 
10 
54 
41 
7.2 
12 
25 
25 
19 
18 
38 
68 
8.0 
7.8 
28 
27 
20 
14 
44 
46 
32 
11 
15 
140 
23 
120 
210 
130 
14 
24 
47 
140 



8.2 
7.4 
7.6 
31 
10 
25 

7.1 
12 
26 
26 



7.1 
7.9 
28 
26 

12 
42 
36 

10 
11 

140 
20 
70 

150 
69 



300 
300 
300 
900 
300 
900 
2.5 
300 
300 
300 
300 
300 
300 
300 
2.5 
300 
300 
300 
300 
300 
300 
300 
300 
300 
300 
300 
120 
300 
120 
120 
120 
180 
900 
300 
900 



0.79 + 0.02 
0.78 + 0.02 
1.96 + 0.04 
0.78 + 0.02 
1.95 + 0.04 



0.90 + 0.02 
2.12 + 0.04 
2.12 + 0.04 

2.12 + 0.04 
0.90 + 0.02 

2.13 + 0.04 



1.26 + 0.03 
2.60 + 0.05 
2.74 + 0.05 

2.64 + 0.05 
1.25 + 0.03 
2.59 + 0.05 
2.73 + 0.05 

2.65 + 0.05 

2.94 + 0.06 
6.03 + 0.12 

2.98 + 0.06 
6.14 + 0.12 

5.99 + 0.12 

6.46 + 0.13 
1.34 + 0.03 
4.10 + 0.08 
0.60 + 0.01 

1.47 + 0.03 



0.81+0.04 3 

0.81 + 0.04 3 

1.84 7 

0.81 + 0.04 3 

1.84 7 



0.98 + 0.12 

2.13 

2.22 

2.18 

0.98 + 0.12 

2.18 



1.27 + 0.02 
2.53 + 0.09 
2.62 + 0.09 
2.57 + 0.09 
1.27 + 0.02 
2.53 + 0.09 
2.62 + 0.09 
2.57 + 0.09 

2.88 + 0.08 
6.44 + 0.24 
2.88 + 0.08 
6.32 + 0.24 
6.44 + 0.24 
6.51 + 0.24 
1.20 + 0.10 
4.08 + 0.18 
0.62 + 0.09 
1.30 



References: 1 iBecker. W hite & Edwar d^ (|l99lh._ 2 ICregorv fc CondonI Jl99lh 3 ICriffith et all (|l995l ). 4 iKiihr et al.l (|l98lh . 5 

I(ll990|). 



.^^. ^ „ ^^ J. 2 |G 

iLaing fc Peacock! (ligSlf Te iWhite fc Beckeil Jl992l ). f lwright fc Otrupcek 



^ Although the 1980 and 1981 observations have a centre frequency of 4885.1 MHz, the weighted mean for the combined dataset is 
still close to 4860.1MHz 



Before combining configurations, we subtracted in the {u, v) 
plane all sources outside a fixed central field. For 0755+37 
at L-band, this procedure failed to remove sidelobes at the 
centre of the field from a bright confusing source close to 
the half-power point of the primary beam. The reason is 
that the VLA primary beam is not azimuthally symmetric, 
so the effective complex gain for a distant source is not the 
same as that at the pointing centre and varies with time 
in a different way. We used the AiPS procedure peelr to 



remove the offending source from the {u,v) data for each 
configuration before combining them. 

Finally, we corrected for variations in core flux density 
and amplitude scal e between observations as described in 
iLainget all (|2006bl ). 

J2000 coordinates are used throughout this paper. If 
positions from archival data were originally in the B1950 
system, then (u,v,w) coordinates were recalculated for J2000 
before imaging. The astrometry for each of the sources was 
set using the A-configuration observations, referenced to a 
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nearby phase calibrator in the usual manner. Thereafter, the 
position of the compact core was held constant during the 
process of array combination. 

The observations of M84 in 1980 and 1981 used an ear- 
lier and less accurate value of the position of the phase cal- 
ibrator B1236+077 (alias J1239+075) than that currently 
given in the VLA calibrator manual. We have updated the 
astrometry to reflect the improved calibrator position. The 
archival L-band observations of M 84 taken in 2000 used a 
pointing centre displaced by ?^1.1 arcmin from the centre of 
the source. 

The C-band data were usually taken in two adjacent 
50-MHz frequency channels, which were imaged together. 
The L-band channels were also imaged together in /, in or- 
der to derive spectral-index images. For all sources except 
0326+39, they were also imaged independently, primarily for 
analysis of linear polarization. 

In order to avoid the well-known problems introduced 
by the conventional CLEAN algorithm for well-resolved, 
diffuse brightness distributions, total-intensity images at 
the higher resolutions were produced using the multi- 
scale CLEAN algorithm a s implemented in the AiPS pack- 
age (jGreisen et al.l l2009ll or, in one case, a maximum- 
entropy alg orithm (ICornwell fc EvansI Il985l , used as de- 
scribed by iLeahv fc Per lev! Il99lh . The standard single- 
resolution CLEAN was found to be adequate for the lowest- 
resolution / images. Stokes Q and U images were CLEANed 
using one or more resolutions (we found few differences be- 
tween single and multiple-resolution clean for these images, 
which have little power on large spatial scales). All of the 
images were corrected for the effects of the antenna primary 
beam. 

In general, the deep 4.9-GHz images have off-source 
noise levels very close to those expected from thermal noise 
in the receivers alone. There are a few faint artefacts on the 
/ images. These are visible as concentric rings around the 
bright cores and, for NGC 193 at 1.35 and 1.6-arcsec res- 
olution, a quadrupolar pattern at the 2a level. These are 
due to errors in cross-calibration of the different array con- 
figurations, which were particularly troublesome due to the 
low declination of this source. The integrations for all of the 
L-band images and for the C-band image of 0326+39 are 
shorter and confusion from sources outside the field of view 
is worse, so noise levels are correspondingly higher. 

Finally, we produced improved / images from self- 
calibrated L- and X-band visibility data for 3C 296 
(|Laing et al.ll2006bl ) using the multi-resolution clean algo- 
rithm. 

As a check on the amplitude calibration and imaging 
of the / images used in spectral-index analysis, we inte- 
grated the flux densities using the AiPS verb tvstat. We 
estimate that our errors are dominated by a residual scale 
error of ~2%. All of the results are in excellent agreement 
with single-dish measurements (Table [3]). 

The configurations, resolutions, deconvolution algo- 
rithms and noise levels for the final images are listed in 
Table \3\ The noise levels were measured before correction 
for the primary beam, and are appropriate for the centre of 
the field. 



2.3 MERLIN observations and reduction 

We also present MERLIN imaging in total intensity only 
for two of the sources. 0206+35 was observed for a total 
time of about 14 hours. The array included the following 
telescopes: Defford, Cambridge, Knockin, Wardle, Darnhall, 
Mk2, Lovell and Tab ley. The observations were carried out 
at 1420 MHz with a bandwidth of 15 MHz, in each of left 
and right circular polarizations. The nearby compact source 
0201+365 was used as the phase calibrator and the flux- 
density scale was determined using 3C 286. The data were 
edited, corrected for elevation-dependent effects and non- 
closing errors and flux-calibrated using the standard MER- 
LIN analysis programs. Imaging and self-calibration were 
again performed using the AIPS package. The off-source im- 
age rms after self-calibration was close to that expected from 
receiver noise alon e. The MERLIN ob servations of 0755+37 
were described bv lBondi et al.l (|2p0d). 

The parameters of both MERLIN images are given in 
Table El 



3 IMAGES 

3.1 General 

Our conventions for Figs [Tl - [T2l and the descriptions in the 
text are as follows. 

(i) Images of total intensity, /, are shown as grey-scales, 
over ranges indicated by the labelled wedges. The units are 
mJybeam"^. 

(ii) We also show grey-scales of intensity gjradient, |V/ |, 
approximated using a Sobel filter (jSobel fc Feldmanlll968l ). 

(iii) We use the notation P = {Q^ -\- U^y^^ for polar- 
ized intensity and p = P/I for the degree of linear polar- 
ization, pi, is the degree of polarization at frequency u (in 
GHz). All values of P hav e been corrected for Ricean bias 
(| Wardle fc KronberdI 19741 ) . Linear polarization is illustrated 
by plots in which vectors with lengths proportional to the 
degree of polarization at 4.9 GHz (^4.9) and directions along 
the apparent magnetic field (Ba) are superposed on false- 
colour images of either / (again with a labelled wedge indi- 
cating the range) or |V/|. A value of p = 1 is indicated by 
the labelled bar. The apparent field direction is xo + vr/2, 
where xo is the E-vector position angle corrected to zero- 
wavelength by fitting to the relation x('^^) = Xo + RMA^ 
for foreground Faraday rotation derived from the images 
in Table [3] (RM is the rotation measure). In some sources, 
we used RM images at lower resolution to correct the po- 
sition angles, as detailed in the captions. This procedure 
is valid if the RM varies smoothly over the low-resolution 
image, and maximises the area over which we can deter- 
mine the direction of the apparent field. Vectors are plot- 
ted where: (a) / > 5(7/ , (b) P > 3a p (the noise levels 
are given in Table [3| and (c) the RM is well-determined. 
The RM images for 0206+35 and M 84 were determined us- 
ing three and fou r frequencies, respectively and are shown 
in Guidet ti et al.1 ([2OII); that for 0755+37 (also from three 
frequencies) will be described by Guidetti et al. (in prepa- 
ration). For NGC 193, images were only available at 4.86 
and 1.365 GHz. T he integrated RM of the source is small 
( — 18 ± 2radm~^; ISimard-Normandin, Kronberg &: ButtonI 
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Figure 1. Images of NGC 193 at 4.05-arcsec resolution, (a) Total intensity at 4.9 GHz. (b) Intensity gradient for the image in panel (a). 
The arrows mark the high brightness gradients which coincide with the edges of the flat-spectrum caps (which are also marked), (c) 
Spectral-index image between 4.9 and 1.4 GHz. Values outside the white contour are lower limits. The caps are indicated, (d) Vectors 
with directions along Ba and magnitudes proportional to P4.9, plotted on a false-colour image of total intensity at 4.9 GHz. 



Il98lh as are the variations of position-angle difference across 
the source. We are therefore confident that a two-frequency 
RM determination is adequate in this case (the correction 
required to derive Ba is in any case very small). 

(iv) Spectral index, a, is defined in the sense S{iy) oc iy~^ 
and we use a^l for the index between frequencies ui and 1^2 
(in GHz). In the false-colour images of spectral index, the 
input / images at the lower frequency are always blanked 
for / < 3a I O In cases where the areas over which emission 
was detected were essentially the same at both frequencies, 
we also blanked the higher- frequency image at 3a i. If sig- 
nificant areas were detected only at the lower frequency, we 
did not blank the higher- frequency image. Instead, we plot 
a single contour which indicates the boundary of the region 
where the source is detected at / > 3ai at the higher fre- 
quency. Outside this contour, the spectral indices are lower 
limits. We have carefully inspected the spectral-index im- 
ages to check for edge effects and zero- level problems. We 
are confident that the values and lower limits in all of the 



^ The spectral-index image for M84 has additional blanking, as 
noted in the caption of Fig. [9{c) 



unblanked regions are reliable except where explicitly stated 
and that the steep spectra seen at the edges of the lobes in 
several sources are real. 

(v) The restoring beam (FWHM) is shown at the bottom 
of each plot. 

(vi) We refer to parts of the sources by the abbreviations 
N, S, E, W (for North, South, East, West) etc. 

(vii) We refer to the main (brighter) and counter (fainter) 
jets. 



3.2 NGC 193 

Fig.[T]^a) shows the total-intensity distribution over NGC 193 
at 4.9 GHz and 4.05-arcsec FWHM resolution. The sym- 
metrical jets appear to broaden rapidly and also bend away 
from their initial straight path as they reach the midpoints 
of symmetric lobes. The lobes both have well-defined lead- 
ing edges that approximate arcs of circles in projection on 
the sky (the W lobe having a larger radius of curvature), 
but they lack hot spots that might mark the termination 
points of the jets. A broad, faint emission bridge fills the 
central region of the source between the lobes and appears 
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Figure 2. High-resolution images of the inner jets of NGC 193. (a) Total intensity at 4.9 GHz with 1.6arcsec FWHM resolution, (b) 
Brightness gradient derived from the image in (a). Sharp steps in brightness ("arcs") are indicated, (c) Spectral index between 4.9 and 
1.4 GHz at 1.6-arcsec resolution. The most prominent of the arcs, which shows a flattening in spectral index, is marked, (d) Vectors with 
directions along Ba and magnitudes proportional to ^4.9, plotted on false-colour images of total intensity at 4.9 GHz with 1.35arcsec 
FWHM resolution. The vectors have been corrected for Faraday rotation using a 4.05-arcsec FWHM resolution RM image, (e) As in 
panel (d), but for the innermost jet regions at 0.45-arcsec resolution. 
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Figure 3. Images of 0206+35. (a) Total intensity at 1.4 GHz with 4.5-arcsec FWHM resolution, (b) Total intensity at 4.9 GHz, 1.2-arcsec 
FWHM. (c) Intensity gradient at 1.4 GHz, 4.5-arcsec FWHM. (d) Intensity gradient at 1.4 GHz, 1.2-arcsec FWHM. The arrows mark the 
high brightness gradients around the inner boundary of the flat-spectrum cap. (e) Spectral index between 4.9 and 1.4 GHz, 4.5-arcsec 
FWHM. (f) Spectral index between 4.9 and 1.4 GHz, 1.2-arcsec FWHM. In panels (e) and (f), values outside the contours are lower 
limits. 



to be wider than the lobes in the N-S direction in the centre 
of the s ource (similar to the "wings" observed in some PR II 
sources: ICheundl2007l . and references therein). Figs[TJb) and 
(c) taken together show that broad "caps" of emission can 
be delineated in both lobes by enhanced intensity gradients 
and by lower-than-average values of a\X. The intensity gra- 
dients are largest at the outer edges of these caps, but there 
are also gradient features within the lobes (marked with ar- 
rows on Fig.llJ)) which coincide with the edges of the flatter- 
spectrum region. The regions where the jets are most promi- 
nent have low spectral indices around 0.6. The spectral in- 
dex at the trailing edges of the caps steepens smoothly to 
a ~ 0.9 where the emission merges with the broader, sym- 
metric lobes. The most diffuse lobe emission has spectral in- 
dices increasing from ?^1 at the edges of the most elongated 
parts of the lobes to ^lA near the centre of the source, a 
spectral-index pattern characteristic of lobed radio sources 
of both FR classes (see Section 14. 2p . There are regions of 



emission with spectral index approaching 2 at the edges of 
the faintest emission on the N and S edges of the sourcefj 
Fig. [Hd) shows that the distribution of the apparent mag- 
netic field direction over the lobes is basically circumferen- 
tial, while the magnetic field in the jets is perpendicular to 
the jets over most of their lengths. The structure of the ap- 
parent magnetic field in both the jets and the lobes appears 
regular, and characteristic of that commonly found in jets 
of FRI sources and in the lobes of both FR classes. 

Fig. [2ja) shows the total-intensity distribution over the 
jets at 4.9 GHz with a resolution of 1.6arcsec FWHM. Both 
jets are fairly straight and similar in overall appearance, ex- 
hibiting rapid lateral expansion just beyond the distance 
from the unresolved nuclear radio source at which the E jet 
is markedly brighter than the W jet. Their edges are well 



^ The spectral index within a few arcsec N and S of the core is 
affected by artefacts in the 4.9-GHz image. 
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Figure 4. Vectors with directions along the apparent magnetic field and lengths proportional to the degree of polarization at 4.9 GHz, 
superimposed on a false-colour image of total intensity across 0206+35 at the same frequen cy. The reso lution is 1.2arcsec FWHM. The 
apparent field directions were derived using a three- frequency rotation-measure fit (jGuidetti et al.ll2011[ ). 



delineated by steep transverse intensity gradients near the 
midpoint of the source (Fig. [2]3). The surface brightnesses 
of both jets decrease smoothly with distance from the nu- 
cleus except at a distance of ?^25arcsec, where there are 
more sudden drops, mostly clearly visible as "arcs" cross- 
ing the jets on the gradient image (indicated on Fig. [2J)). 
The overall spectral index of the jet-dominated emission ap- 
pears to steepen with distance from the nucleus, but this is 
almost certainly the result of superposition of dimming jets 
on steeper-spectrum lobe emission. The prominent arc in the 
E jet is associated with a slight flattening in the spectrum 
(Fig.[2t). 

Fig. EJd) shows the intensity and apparent magnetic 
field distributions over the inner ?^45-arcsec regions of both 
jets at 1.35-arcsec resolution. The magnetic field organiza- 
tion over both jets is quite regular, with the field closest to 
the jet axis predominantly orientated perpendicular to the 
axis. The apparent field at the edges of the inner jets is par- 
allel to the rapidly-expanding outer isophotes. Farther from 
the nucleus, the edge field directions in both jets converge 
towards the axis to form almost circular patterns. 

Fig.[2je) shows the 4.9-GHz total intensity and apparent 
magnetic-field distributions over the inner 15 arcsec (?^4.5 
kpc) of the jets at 0.45-arcsec resolution. Both jets exhibit 
faint inner regions in the first ~2 arcsec from the nucleus 
before they brighten and subsequently flare. There is bright, 
non-axisymmetric knot structure in the fl rst ^4 arcsec of th e 
E jet, downstream of the flaring point (|Laing et al.lll999h , 
where both jets brighten abruptly. 

3.3 0206+35 

Fig. [3] shows the total-intensity, brightness gradient, and 
spectral-index distributions over the whole of 0206+35 at 



two resolutions. Fig. (S^a) at 1.4 GHz, 4.5-arcsec FWHM 
resolution, shows that the large scale structure consists of 
two lobes, overlapping and circular in cross-section, with 
well-defined outer edges to the NW and SE of the source, 
superimposed on fainter diffuse emission to the N and S. 
Fig. [3jb), at 4.9 GHz and 1.2-arcsec resolution, shows the 
source in more detail but with reduced sensitivity to the 
largest-scale emission. At this resolution, both lobes show 
sharp outer boundaries. The roughly circular edge of the 
NW lobe protrudes beyond the diffuse emission, whereas 
the corresponding feature in the SE is well inside the outer 
boundary of the source and is most obvious in the E of the 
lobe, close to the termination of the jet. If the orientation 
of ^ ~ 40° determined for the inner jets (Laing & Bridle, 
in preparation) also applies to the lobes, then they are pre- 
sumably ellipsoidal with an axial ratio ?^1.6. Fig. [S^b) also 
shows some internal structure in both jets. The NW jet has 
the brighter base, and both bends and brightens as it enters 
its lobe, after which its path meanders. The SE counter-jet 
appears to expand more rapidly initially, then also meanders 
as it enters its lobe. 

Figs [3{c) and (d) show the 1.4-GHz intensity gradient 
images of 0206+35 at 4.5-arcsec and 1.2-arcsec resolution, 
respectively. The edges of the jets are clearly marked by 
enhanced intensity gradients at both resolutions, while sig- 
nificant internal structure is also apparent in the lobes. Both 
lobes exhibit strong brightness gradients at their outer edges 
in these displays, corresponding to the sharp boundaries 
noted earlier. There is a particularly striking correlation be- 
tween the main features of these intensity-gradient images 
and of the two 1.4 to 4.9-GHz spectral- index images shown 
as Figs. [3{e) and (f). The emission inside the brightest in- 
tensity gradients around the jet has a much lower spectral 
index, typically <0.65, than the ?^0.7 to 1.0 spectral index 
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Figure 5. (a) Grey-scale of the 4.9-GHz total-intensity distribu- 
tion over the jets in 0206+35 at 0.35-arcsec FWHM resolution. 
The grey-scale range is - 2.5 mJy beam"-*^. (b) Vectors with 
lengths proportional to the degree of polarization at 4.9 GHz and 
directions along the apparent magnetic field, superimposed on a 
false-colour display of the total intensity at 4.9 GHz. The resolu- 
tion is 0.35 arcsec FWHM. The vector directions are derived from 
3-frequency RM fits at 1.2-arcsec resolution, (c) Grey-scale of the 
1.6-GHz total intensity distribution over the inner 5 arcsec of the 
NW jet and the unresolved nuclear source in 0206+35 at 0.16- 
arcsec FWHM resolution, from MERLIN. The grey-scale range is 
0-5 mJy beam" ^ . 



that is prevalent over the rest of the two spherical lobes. The 
diffuse emission outside the lobes has spectral indices rang- 
ing from ?^1.05 to >2, generally increasing with distance 
from the lobes towards the outer edge of the source. Also 
notable are the "fans" of lower-spectral-index emission that 
can be traced from the ends of the jets to the regions at the 
edges of both lobes that show the most pronounced bright- 
ness gradients. This is particularly striking in the NW lobe, 
where a cap of lower-spectral-index emission is bounded by 
the high brightness gradients marked by arrows on Fig. [3jd) 
and by the edge of the lobe. This suggests that the jet out- 
flow has reached the end of the lobe in a less-collimated, but 
still identifiable, form. In the SE lobe, the jet bends to the N 
before appearing to impact the edge of the lobe at another 
enhanced brightness gradient, marked in Fig. [3{d). 

Fig. m shows that the magnetic field configuration in 
both lobes is well ordered and basically circumferential, 
while the magnetic field in the jets is predominantly perpen- 
dicular to their axis near the centre lines of the jets, with 
evidence for parallel field at the jet edges. The southern edge 
of the source is strongly polarized with field tangential to the 
boundary. 

At 0.35-arcsec FHWM resolution (Fig. [5ji) the lobe 
emission is substantially resolved out so the images are dom- 
inated by the jets. The bright base of the main (NW) jet is 
clearly centre-brightened, while the corresponding segment 
of the counter-jet appears centre- darkened. While at first 
sight the main jet appears to expand more slowly than the 
counter-jet, the geometry of the centre- darkened segment of 
the counter-jet is very similar to that of the main jet over 
the first ^10 arcsec. This suggests a two-component view of 
these jets wherein a narrow inner structure, brighter to the 
NW and fainter to the SE, is seen superposed on a broader 
expanding structure that is slightly brighter to the SE than 
to the NW. We will interpret this elsewhere (Laing & Bri- 
dle, in preparation) as evidence for a symmetrical relativist ic 
outflow surrounded by a mildly relativistic backflow in this 
source. 

The magnetic field is clearly perpendicular to the jet 
axis over most of the length of both jets (Fig. [SJ)), but the 
first few arcsec of the main jet, where it is brightest, have the 
magnetic field parallel to the jet axis. There is also evidence 
for oblique, or parallel, magnetic field at the edges of both 
jets. 

Fig. E^c) shows the bright, narrow base of the NW 
jet from MERLIN data at 1.6 GHz. This higher-resolution 
(0.16 arcsec FWHM) image unambiguously identifies the 
flaring point in the jet, 0.7 arcsec from the nucleus, where 
it brightens abruptly. This is an important fiducial distance 
for our modelling. The image also shows the start of rapid 
expansion downstream of the flaring point. 



3.4 0755+37 

Fig. [6] shows the total-intensity, |V/| and a distributions 
over all of 0755+37 at two resolutions and the polarization at 
low resolution. Fig. [6la) at 1.4 GHz, 4.0-arcsec FWHM res- 
olution, shows that the large-scale structure consists of two 
lobes, again roughly circular in projection, with well-deflned 
but not particularly sharp outer edges to the W and E, plus 
fainter diffuse emission to the N and S. The E lobe has a se- 
ries of narrow ridges and brightness steps, all roughly arcs of 
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Figure 6. Images of the whole of 0755+37 at resolutions of 4.0arcsec FWHM (panels a, c, e and g) and 1.3arcsec FWHM (panels b, d 
and f). (a) Total intensity at 1.4 GHz. (b) Total intensity at 4.9 GHz. (c) Intensity gradient at 1.4 GHz. (d) Intensity gradient at 4.9 GHz. 
(e) and (f) spectral index between 4.9 and 1.4 GHz. (g) vectors with lengths proportional to p4.9 and directions along the apparent 
magnetic field from a three- frequency rotation-measure fit (Guidetti et al., in preparation), superimposed on the intensity gradient at 
4.9 GHz. 
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circles in projection, in the region where the brighter jet ap- 
pears to terminate. They are recessed from the E boundary 
of this lobe and some may be the edges of thin shells. The 
structure of the W lobe is unusual, containing some arc-like 
features and other structure suggestive of a rapidly decolli- 
m atingi; counter - jet W of the nucleus, as previously described 
bv lBondi et al.l (|2000l ) . with a "hole", or deficit of emission in 
the region where the counter-jet might be expected to ter- 
minate. Fig.ElJb) at 4.9 GHz, 1.3-arcsec FWHM resolution, 
is insensitive to the largest scales of emission to the N and 
S of the main source, but clearly shows internal structure in 
the E lobe, including the concentric semicircular ridges at 
the end of the jet (labelled on the figure). All of the sub- 
structure in the W lobe appears to be resolving out, though 
vestiges of the ridge apparent at lower resolution in Fig. [6ja) 
remain. 

Figs. [6jc) and (d) show the intensity gradients over the 
whole source at 1.4 GHz, 4.0arcsec FWHM resolution and 
4.9 GHz, l.Sarcsec FWHM resolution, respectively. These 
figures emphasise the strong differences between the internal 
structures of the lobes: multiple recessed ridges with signifi- 
cant brightness gradients in the E lobe, but much smoother 
structure in the W lobe away from the jet. 

Figs. EJe) and (f) clearly show three distinct spectral- 
index regions on each sides of the source, as follows. 

(i) a ~ 0.6 at the bases of both jets, in a broad cap in the 
NW part of the W lobe, and all along the region delineated 
by the strongest brightness gradients in the E jet. 

(ii) a ^ 0.8 over most of the rest of both lobes, including 
the ridge extending Northward from the nucleus. 

(iii) There is also steeper-spectrum diffuse emission with 
a increasing from ?^ 1 in the central part of the source to 
~1.5 at the N and S edges. 

The spectral-index image suggests that the counter-jet flow 
persists as far as the ridge of emission in the W lobe marked 
on Figs.[6ja) and (b), despite the lack of evidence for this in 
total intensity. 

Fig.[6{g) shows that the apparent magnetic field in both 
lobes is exceptionally well organised, and mainly tangential 
to the lobe boundaries. The degree of linear polarization is 
p^O.6 over much of both lobes, consistent with the high 
degree of organization evident from the vectors. Note the 
excellent alignment between the field vectors and the ridges 
of high brightness gradient on both sides of the source. Fig. [71 
confirms the exceptional degree of ordering of the magnetic 
field in both lobes at 1.3-arcsec FHWM resolution. 

Fig. [8] shows the jet bases on larger scales and at higher 
resolution. Fig. [8ja) is optimised to emphasise fine-scale 
structure in the jets at 1.3-arcsec resolution. It shows that 
the E jet has the brighter base but becomes limb-brightened 
at the position indicated on Fig.[8ja) at about 18 arcsec from 
the nucleus. At this resolution, the W counter-jet contains a 
centre- darkened structure that expands at about the same 
rate as the brighter E jet, embedded in a much broader, 
and more rapidly-expanding cone of emission with at least 
two curved arcs (Fig. [8^). The centre of the counter-jet is 
crossed by a prominent, straight bar of emission (labelled 
as such on Fig. ^,) at ?^15 arcsec from the core. The steep- 
est brightness gradients at the edges of both jets are very 
similar in form within ?^ 15 arcsec of the nucleus so that, as 
in 0206+35, the mner geometry of the counter-jet structure 



appears to mimic that of the brighter jet on the other side 
of the source (Fig. (Sj)). 

Fig. [8{b) also shows detail of the magnetic field organi- 
sation at the bases of both jets at 1.3-arcsec resolution. The 
bright base of the E jet has the magnetic field roughly par- 
allel to the jet axis, but there is a rapid transition, with the 
field becoming perpendicular to the axis as the jet expands. 
On the counter-jet side the field is also transverse. There 
is little evidence for any perturbation of the magnetic field 
structure at the edges of the jets except at the N edge of 
the counter-jet where the magnetic field becomes parallel to 
the steepest brightness gradient once the jet widens signifi- 
cantly. The very high degree of polarization in the surround- 
ing diffuse emission makes it difficult to disentangle the true 
polarization of the jets where their emission is weak, e.g. at 
their edges. 

Figs [8{a) and (b) also show a filament of faint emis- 
sion which extends for about 30 arcsec Northward from the 
vicinity of the nuclear source, roughly perpendicular to the 
jets. It has a very high degree of linear polarization, with 
an apparent magnetic field parallel to its length. Fig. [7] sug- 
gests that this highly polarized filament may be part of a 
larger region of enhanced polarization that delineates the 
inner boundary of the W lobe. 

Fig. [8{c) shows the total-intensity structures at the 
bases of the jets at 0.4 arcsec resolution. The main jet is 
clearly centre-brightened whereas the counter-jet is not and 
the brighter edges of the counter-jet lie mostly outside the 
region that would be delineated by reflecting the main jet 
across the nucleus. As in 0206+35 (Fig.[5K) there is a narrow 
collimated structure within which the main jet is systemati- 
cally brighter than the counter-jet apparently superposed on 
a broader structure which is brighter around the counter-jet 
than around the main jet. As for 0206+35, we will show 
elsewhere (Laing & Bridle, in preparation) that this struc- 
ture can be modelled as a symmetrical relativistic outflow 
surrounded by modestly relativistic backflow. 

The polarization image in Fig. [8{d) shows the extent 
of the region at the bright base of the main jet in which 
the magnetic field at the edges is parallel to the expanding 
outer isophotes, whereas the on-axis field is oblique (see also 
Fig. [SJ)). Finally, a 1.7-GHz MERLIN image of the main jet 
base (Fig. [8^) shows the position of the flaring point and the 
details of the initial expansion. 



3.5 M84 

M 84 is of particular interest for two reasons: it has a much 
lower radio luminosity than the other sources we have stud- 
ied and it shows v ery clear evidence f or interaction with the 
surrounding IGM (JFinoguenov et al.| [2008). 

Fig. (Ha) shows the total-intensity distribution over 
M84 at 4.9 GHz with 1.65-arcsec FHWM resolutior0. Both 
jets of this small (overall extent ~12 kpc), low- luminosity 
radio source expand rapidly and deflect within about 1 ar- 
cmin. They are surrounded by diffuse emission (at least in 
projection) everywhere except perhaps within a few arcsec 
of the nucleus. The initially brighter N jet can be traced 



^ Lower-resolution (FWH M ?:i4 arcsec) radio observations, shown 
bv lLaing fc Bridla (|l987l ), are not reproduced here. 
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Figure 7. Vectors with lengths proportional to the degree of polarization at 4.9 GHz and directions along the apparent magnetic field, 
superimposed on a false-colour display of the total intensity over 0755+37 at 4.9 GHz. The resolution is 1.3 arcsec FWHM and the colour- 
scale range is - 2.5 mJy beam"-*^. The vector directions are derived from three-frequency RM fits at 4.0-arcsec resolution (Guidetti et 
al., in preparation). 



as far as the edge of its lobe, where it bends through ?^90° 
in projection and decollimates on impact. The bending is 
accompanied by strong brightness gradients (Fig. [SId). In 
contrast, and uniquely amongst the sources in this paper, 
the S jet (initially fainter and misaligned with the nucleus) 
appears to terminate within its lobe and to feed a bubble- 
like structure with significant internal brightness gradients 
and filaments. The bubble is contained within a smoother, 
more elongated structure, at least in projection. The spec- 
tral index between 1.4 and 4.9 GHz (Fig. [St) is constant with 
a ?^ 0.6 over the jets, within the Southern bubble and over 
most of the N lobe. The only regions of significantly steeper 
spectral index that we have detected (marked by arrows on 
Fig. [91:) are on both sides of the south jet base and to the 
west of the north jet. The current observations are too noisy 
to determine the spectral index in the low- surface- brightness 
emission outside the southern bubble. Fig. [H^d) shows the 
apparent magnetic field structure over the whole source at 
1.65-arcsec FHWM resolution. The magnetic field in the S 
lobe is broadly circumferential and appears well-aligned with 
the peak in brightness gradient at the edge of the bubble. 
There is a sudden increase in the degree of polarization at 
the edge of the bubble, suggesting a discontinuity in the 
field structure at that location. A configuration in which the 
field is confined to ellipsoidal shells but is otherwise ran- 



dom (model A of lLaindll980l l gives qualitatively the correct 
polarization distribution, but the predicted variation of p 
across the lobe is smoother than we observe. The magnetic 
field in the N lobe is predominantly perpendicular to the 
presumed path of the jet along its mid-line. 



Fig. llOf a) shows the 4.9-GHz total-intensity distribu- 
tion over the whole source at 0.4-arcsec FWHM resolution. 
This highlights the filamentary structure in the N lobe, the 
S bubble and a thin rim of emission around the S lobe. The 
intensity and gradient images of the jets at this resolution 
(Fig. llOb and c) emphasise the edges of both jets and the 
curved arcs in the north. The former also shows a curious 
thin feature (labelled A) joining the S jet and the edge of its 
lobe. The misalignment (non-collinearity) of the axes of the 
N and S jets beyond a few arcsec (a few hundred parsecs) 
from the nucleus and the initially knotty structure of the N 
jet, can be seen on a larger scale in Fig. [TOT d), which also 
shows faint emission close to the nucleus on both sides. 

Fig. [TOT e) and (f) show that, although the apparent 
magnetic field in the jets is locally well-organised, there are 
significant regions where the field is oblique to the jet axis. In 
the outer parts of both jets the field appears to be predom- 
inantly perpendicular to the jet axis, but the jet emission 
also becomes blended with that from the lobes. 

M84 may be an intermediate case between lobed and 
tailed sources, showing some characteristics of each class. 
The N jet terminates in a sharp bend at the outer edge of its 
lobe, as often seen in lobed sources (eg. 3C 296, Section [X6)) , 
but there are hints of a nascent tail struct ure in the NE. This 
is supported by Chandra imaging of M 84 (Finoguenov e t al.l 
2008), which suggests that the NE lobe is breaking out of 
the surrounding hot plasma. In contrast, the S jet termi- 
nates well within its bubble-like lobe. The oscillation of the 
S jet prior to its eventual disruption is very reminiscent of 
that of the jets within the S spurs of the tailed sources 30 31 
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Figure 8. High-resolution images of the inner jets of 0755+37. (a) Total intensity at 4.9 GHz with 1.3-arcsec FWHM resolution, plotted 
with a compressed grey-scale range to emphasise fine-scale structure in and around the jets, (b) Vectors with lengths proportional to ^4.9 
and directions along the apparent magnetic field superimposed on a false-colour image of intensity gradient at 4.9 GHz. The resolution 
is l.Sarcsec FWHM. (c) Total intensity at 4.9 GHz, 0.4-arcsec FWHM. (d) Main jet base at 4.9 GHz with 0.4-arcsec FWHM resolution. 
Ba vectors with lengths proportional to ^4.9 are superposed on a false-colour plot of total intensity, (e) MERLIN image of the main 
jet base at 1.7 GHz with 0.16-arcsec FWHM resolution (Bondi et al.ll2000l ). Corrections for Faraday rotation in panels (b) and (d) were 
made using a three- frequency RM fit at 4.05-arcsec resolution (Guidetti et al., in preparation). 
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Figure 9. Images of M84 at 1.65-arcsec FWHM resolution, (a) 4.9-GHz total intensity, (b) 4.9-GHz intensity gradient, (c) Spectral index 
between 1.4 and 4.9 GHz, plotted only where its rms error is <0.1. The vertical "streaks" are artefacts. Arrows mark the areas when the 
spectral index is significantly steeper than the typical value oi a = 0.6. (d) Ba vectors with lengths proportional to P4.9, superposed on a 
false-colour plot of in tensity gradient. Corrections for Faraday rotation were made using a 4-frequency RM image at 4.5-arcsec resolution 
(|Guidetti et al.ll201ll ). All panels show identical areas. 
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Figure 10. 4.9-GHz images of M84 at 0.4-arcsec FWHM resolution, (a) Total intensity for the whole source, (b) Total intensity for the 
jets, (c) Intensity gradient for the jets. Three prominent "arcs" in the N jet are labelled on panels (b) and (c). (d) Total intensity for the 
inner jets, showing the abrupt bend in the counter-jet. (e) Ba vectors with lengths proportional to P4.9, superimposed on a false-colour 
image of intensity gradient for the N jet. (f) As (e), but for the S jet. 
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and 3C 449 (JLaing et al.ll2008l : lKatz-Stone fc Rudnickll997l l. 
The spectral gradients (such as they are) are more character- 
istic of lobed sources, with no hint of steepening outwards. 
The constancy of the spectral index across the radio struc- 
ture is not surprising: given the small (^12-kpc) size of the 
source, synchrotron losses are unlikely to have had enough 
time to steepen the spectrum at GHz frequencies even in the 
more extended regions. 



3.6 3C296 

Fig.llllshows the 1.4-GHz total intensity and brightness gra- 
dient together with the 1.4 to 4.9-GHz spectral index distri- 
butions over 3C296 at 5.5-arcsec resolution. T he intensity 
data are essentially those in Figs. 1(a) and 4(a) of Lain g et al.l 
([2006b) with an improved deconvolution, but the grey-scale 
range in Fig. [TIT a) is chosen to show the jets more clearly 
where they appear to enter the lobes. The corresponding in- 
tensity gradient is shown in Fig. [TTJb). Lower limits to the 
spectral indices have been plotted outside the white inten- 
sity contour in Fig. \TT\ c) to provide a better representation 
of the large-scale spectral gradients at the edges of the ra- 
dio source. There is clear evidence that the flatter-spectrum 
{a ^ 0.5 to 0.65) jets propagate to the edges of both lobes, 
where they deflect and eventually blend with more extended 
steeper-spectrum emission whose spectral index a ^ 1. The 
NE jet forms a cap of flat-spectrum emission with a semi- 
circular outer boundary, again marked by sharp brightness 
gradients. The flow (as traced by its flatter spectrum) then 
turns through ^ 140° in projection, crosses the entire lobe 
and impacts on the boundary at the position marked on 
Fig. [TTr bV The SW jet, on the other hand, does not form a 
cap, but appears to make an oblique impact on the wall of 
the lobe before turning through almost 180° in projection 
back towards the nucleus. The spectral index of the more 
extended emission increases further towards the centre of 
the source and towards its outer edges, where a ~ 2. 



3.7 0326+39 

Fig.ll2lshows the distributions of the 1.4-GHz total intensity 
and 1.4 to 4.9-GHz spectral index over 0326+39 at 18-arcsec 
FWHM resolution. 

As in the other sources studied here, the lobes of 
0326+39 appear to surround the jets in projection. Even at 
this relatively low resolution, the jets are clearly traceable 
to the outer parts of the lobes in both total intensity, where 
they appear to twist and deflect close to the outer edges of 
the lobes, and spectral index, where they can be traced as 
regions of significantly lower index. The W jet exhibits a 
particularly strong kink towards the S about 2 arcmin (^ 
60 kpc) from the nucleus; this kink is clearly replicated in 
the spectral index distribution. 

The extended emission of both lobes also shows a well- 
defined spectral index gradient, increasing towards the nu- 
cleus from a ^ 1 near the broad caps that appear to be dom- 
inated by the outer jet emission to a significantly steeper 
spectrum with a ^ 2 near the centre of the source. The 
spectral index also appears to increase towards a ^ 2 in the 
outer part of the faint Southward extension of the E lobe. 
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Figure 11. (a) Grey-scale of the 1.4-GHz total intensity over 
3C 296 at 5.5-arcsec FWHM resolution, (b) Intensity gradient im- 
age at the same frequency and resolution as panel (a), (c) False- 
colour plot of the 1.4 to 4.9 GHz spectral index distribution at 
5.5-arcsec FWHM resolution; data plotted outside the white con- 
tour are lower limits. 
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Figure 12. (a) Grey-scale of the 1.4-GHz total intensity over 
0326+39 at 18-arcsec FWHM resolution, (b) False-colour plot of 
the spectral index (q;|"^) distribution at the same resolution. 



4 DISCUSSION 

4.1 Initial jet propagation 

There appear to be few, if any, morphological differences 
between the jet base regions in lobed and tailed PR I sources, 
whose common properties include the following. 

(i) The initial rapid expansion (flaring) and recollimation 
of the jets is essentially identical in both types of PR I source. 

(ii) Jet bases usually show significant side-to-side asym- 
metries, although a few very symmetrical examples of each 
type are known. 

(iii) With the exception of these symmetrical cases, there 
are further common properties, as follows. 

(a) One jet in each source exhibits a bright region at 
its base, often with with non-axisymmetric knots and a 
predominantly longitudinal magnetic field. 

(b) Jet brightness and polarization asymmetries are 
correlated: the apparent magnetic field on- axis in the 
brighter jets is initially longitudinal, but switches to trans- 
verse at larger distances; that in the fainter jets is always 
transverse. 

(c) The jet /counter- jet ratio falls with increasing dis- 
tance from the nucleus. 

These regions at t he bases of PR I jets are also those in 
which our models ([Laing &: Bridld l2002al : ICanvin k, Laind 



l2QQ4l : ICanvin et aDl2QQ5l : iLaing et al.ll2006bh show that the 
jets decelerate from relativist ic to sub-relativist ic velocities. 
The development of a transverse velocity gradient across the 
jets implies that they decelerate primarily by boundary-layer 
entrainment of the external medium. We suggest that this 
entrainment occurs primarily in the dense, kpc-scale coro- 
nae of hot plasma that surround the nuclei of twin-jet radio 
galaxies (probably with an additional contribution from stel- 
lar mass loss) and that the entrainment effectively turns off 
on large scales, as it must to avoid further decollimation. 
Chandra observations h ave revealed co ronae of this type in 
0206+35, 0755+37 (Wo rrall et al.ll200lh. 0326+39 (Hardcas - 
tle, private communicat ion), M84 (IFinoguenov et a l.**2008h 
and 30 296 ( Hardcastle et al.1 l2005l : lOroston et al.l .2008:): 
NGO 193 may have a similar component (|Giacintucci et al.l 
|201l|) . The coronae for which data are available have central 
electron densities of 10^ to 7 x 10^ m~^, central pressures of 
3 X 10"^^ to 4 X 10"^° Nm"^ and core radii of 0.3 to 2kpc. 
It seems likely that lobe plasma is excluded from the imme- 
diate vicinity of the nucleus by the high-pressure coronae, 
so the jets in both types of source initially propagate in es- 
sentially identical environments, unshielded from the IGM. 
Further evolution of the jets will depend on their surround- 
ings: if they propagate through low-density lobe material, 
then entrainment will effectively cease and they will recol- 
limate to become almost cylindrical flows, as seen in the 
sources described here. Entrainment rates are also likely to 
be small in jets without surrounding lobes provided that the 
external density is low. For example, there is no sign of any 
lobe surrounding the inner jets of NGO 315, yet it has a ver y 
small opening angle after recollimation (Oanvin et al. 2005). 
This argues for a negligible entrainment rate at distances 
^ 35 kpc, and t he external density d oes indeed fall rapidly 
on these scales ([Oroston et al.ll2008l V In contrast, we have 
argued that entrainment continues at a lower, but still signif- 
icant rate after recollimation in 30 31, whose inner jets also 
have no surrounding lobes. In this source, the opening angle 
after recollimation is larger, our models indicate continuing 
deceleration and there is a hot-gas compon ent with a large 
core radius in addition to the inner corona ([Laing fc Bridld 
l2002b|). 

One feature of jet propagation appears to be unique 
to lobed sources, however. Our high-resolution data for 
0206+35 and 0755+37 show that the apparent difference 
in opening angle between the main and counter-jets seen 
at lower resolution in these sources is a manifestation of a 
two- component jet structure. In both sources, the main jet 
and the counter-jet appear to contain both narrow (well- 
collimated) and broader features on both sides of the nu- 
cleus, but the better collimated parts of the main jet are 
centre-brightened while those in the counter-jet are centre- 
darkened. The broader features at the edges of the counter- 
jets are also slightly brighter than those of the main jets. 
What appears to be poorer collimation of the counter-jet 
at low resolution is now seen as a narrow centre-brightened 
jet opposite a similarly narrow centre- darkened counter-jet, 
surrounded by broader emission which is slightly brighter 
around the counter-jet than around the main jet. We will 
explore explanations for this "two-component" jet structure 
in terms of relativistic outflow in the well collimated com- 
ponent surrounded by mildly relativistic backflow in the 
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broader component in a later paper (Laing & Bridle, in 
preparation) . 

4.2 Jet termination and lobe structure 

With the partial exception of M84, to which we return 
at the end of this section, the sources described in this 
paper have lobes similar t o those in FRII sources (e.g. 
Alexander and Leahvl[l987l: ICarilh et ai]ll99ll : iKharb et all 



20081 ) . They exhibit sharp brightness gradients at their outer 



edges, spectral indices that steepen towards the nucleus from 
the outer lobes and towards the outer extremities of any 
lobe, off-axis wings of diffuse emission near their centres, 
and generally circumferential magnetic fields. The only ob- 
vious difference from the lobes of FR II sources is the evident 
lack of hot spots at the ends of the FRI jets. Similar spec- 
tral gradients occur in a larg er sample of lobed FRI sources 
observed at lower resolution (jParma et al.ll 19991 ) . Where the 
jets dominate the total intensity, they all have similar spec- 
tral indices in the range 0.5 to 0.7. Even when the jets are not 
obviously dominant in intensity alone, the observed spectral 
index distributions can trace plausible paths for jets towards 
the edges of the lobes. This spectral signature implies that 
steeper-spectrum lobe material has been displaced by flatter- 
spectrum jet material along an extended pathway through 
the lobes. Note that the spatial resolution of our data rela- 
tive to the overall source size is higher than for most pub- 
Ushed multifrequency imaging of sources in either FR class. 
Together with the ability to trace the jet flow via its flatter 
spectrum and the use of intensity gradient images to indi- 
cate enhanced compression, our data allow us to identify a 
number of new types of structures in the lobes. 

The jets often terminate in what we have called "caps", 
with the following properties. 

(i) They typically occur at the ends of lobes (one exam- 
ple, 0755+37SE, appears recessed, perhaps as a result of 
projection) and are clearly associated with jet termination. 

(ii) They are bounded at their leading edges by smooth 
outer isophotes (approximated by segments of circles) with 
sharp intensity gradients. 

(iii) They also have inner boundaries, again marked by 
high intensity gradients. 

(iv) Their emission has a flat spectrum, close to that typ- 
ical of jets {a ^ 0.6) after accounting for contamination by 
steeper-spectrum diffuse emission. 

(v) Four out of five examples are fairly symmetrical with 
respect to the local jet axis. The exception is 3C 296NE. 

There are five clear cases of such "caps" out of the ten 
FRI lobes we have studied at high resolution: NGC 193E 
and W, 0206+35NW, 0755+37SE and 3C 296NE, and 
0755+37NW may well be similar. 

Other jet terminations show some, but not all, of the 
same features. In particular, several show enhanced bright- 
ness gradients, but without obvious inner boundaries. In 
0206-h35SE, the jet bends away from its initial direction 
and creates a sharp brightness gradient where it impacts 
on the side of the lobe; the associated emission again has 
a relatively flat spectrum. Both lobes of 0326+39 probably 
have similar structures, but the available resolution is not 
yet high enough to be sure. In 3C 296SW and M84N, the 
jets remain straight until they make oblique impacts on the 



lobe walls at locations marked by high intensity gradients 
after which they bend abruptly. The former case also shows 
flatter-spectrum emission at the impact point. 0755+37NW 
has a weak ring of emission surrounding a flatter-spectrum 
region, but this is contained entirely within the outline of the 
lobe. It is possible that some of these structures (especially 
the last) are caps seen from an unfavourable angle. 

We see little convincing evidence for enhancements in 
I V/| crossing the jets close to their termination points, such 
as might be expected from strong shocks. This implies that 
the flow is internally sub- or transonic. Our estimates of 
the on-axis flow velocities after the initial rapid deceler- 
ation range from ^O.lc to ~0.6c, compared with an in- 



ternal sound speed of 3 



-1/2, 



0.58c for an ultrarelativis- 



tic plasma. Unless there is significant deceleration on scales 
larger than we model, the implication is that the jets must 
be very light and energetically dominated by relativistic par- 
ticles and magnetic field. 

The lobes in this class of source are very different from 
the subsonic, buoyant plumes which are thought to form the 
outer structures of large FRI sources like 3C 31 ([Laing et al.l 
2008). The picture that emerges for jet termination in lobed 
FRI sources is that the flow can be traced at least as far as 
the end of the lobe via its flatter spectrum. Where it impacts 
on the lobe surface, a high-pressure region (the cap) can be 
created. The smooth shape of the outer isophotes (compared 
with the more ragged outline of the lobes) suggests that the 
forward expansion is at least mildly supersonic with respect 
to the external medium. Jet material flows through the cap 
and back into the lobe, eventually mixing with pre-existing 
lobe plasma. The flow pattern is sometimes consistent with 
axisymmetry (or at least appears so in projection) but can 
bend by large angles without completely losing its colli- 
mation. The clearest example is 3C 296NE, where the flow 
bends by ?^140° in the plane of the sky and can be traced 
as far as the trailing edge of the lobe, where its impact is 
marked by a sharp brightness gradient. An alternative to the 
formation of a cap appears to be an oblique collision with 
the boundary of the lobe. In at least one case, 0206+35SE, 
the jet deflects from its original straight path before hitting 
the edge of the lobe. This raises the possibility that the im- 
pact point moves around the surface of the lobe, extending 
it in different directions at different times and lowering the 
average advance speed of the lobe compared with the instan- 
taneous speed of the impact point, as in the "dentist's drill" 
model of FRII sources (iScheuerlll982r ). 



As noted in Section [3.51 M 84 appears to be an interme- 
diate case, in which only one of the jets appears to impact 
on its lobe boundary, but no tails have (yet) developed. The 
S lobe of M 84 is morphologically very similar to the spurs 
in 3C 31, albeit on a much smaller linear scale. This adds to 
the developing picture of the transition between fast, well- 
collimated jets and slow plumes in tailed FRI sources, which 
is often a two- stage process. The jets first enter bubbles, 
within which they disrupt, often thrashing around (as in 
M84S) before disintegrating completely. The tails are then 
formed by escape of material from the bubbles along the 
direction of steepest pressure gradient rather than directly 
from the jet flow. A si milar morphology is often fo und in 
wide-angle tail sources flHardcastle &: Sakellioull2004l V 
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5 SUMMARY 

We have presented deep, high-resolution, multi- 
configuration VLA imaging of four FRI radio sources: 
NGC193, 0206+35, 0755+37 and M84, together with 
lower- resolution observations of 0326+39 and a reanalysis 
of our published images of 3C 296. These sources are 
all examples of "lobed" FRI radio galaxies. Our results, 
displayed as images of total intensity, brightness gradient, 
degree of polarization, apparent magnetic- field direction 
and spectral index, show common features, as follows. 

(i) All of the sources have twin radio jets, with side-to- 
side brightness ratios decreasing with distance from the nu- 
cleus in a manner qualitatively consistent with relativistic, 
decelerating flow. 

(ii) The brightness and polarization distributions of the 
inner jets are very like those in tailed radio sources, indi- 
cating similar deceleration physics. We suggest that the jets 
in both classes of source propagate unshielded from the sur- 
rounding IGM within dense, kpc-scale coronae, leading to 
efficient boundary- layer entrainment. 

(iii) Farther from the nucleus, the jets in both classes of 
source recollimate. This implies that the entrainment rate is 
low whether or not they are surrounded by lobe plasma. 

(iv) 0206+35 and 0755+37 show evidence for a two- 
component jet structure in which a centre-brightened main 
jet and centre- darkened counter- jet are surrounded by 
broader features that are somewhat brighter on the counter- 
jet side, suggesting that a central relativistic outflow is sur- 
rounded by a slower, but still mildly relativistic backflow. 

(v) In all but one case (M84S), the jets propagate to the 
ends of their lobes. Continuing, but less well collimated flow 
can often be traced in spectral index or brightness gradient 
images. 

(vi) Five or six of the ten jets we have studied at high 
resolution terminate at the ends of their lobes in features 
we call "caps" with smooth outer isophotes, sharp inner and 
outer intensity gradients and relatively flat spectra. 

(vii) An additional three out of ten jet terminations are 
best described as oblique collisions of jets with the outer lobe 
walls: they also show enhanced outer intensity gradients and 
flat spectra and may be caps seen at unfavourable angles. 

(viii) The lobes resemble those in FRII sources, with 
sharp outer brightness gradients, spectral indices which 
steepen towards the nucleus and circumferential apparent 
magnetic flelds. 

(ix) There is little evidence for features in the jet bright- 
ness distributions which can be identifled as strong shocks, 
either at recollimation or where they terminat4j. This im- 
plies that the flow is internally sub- or trans-sonic on large 
scales. 

We will present quantitative modelling of the inner jets in 
later papers. 
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